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Largely due to its role in nucleation and crystal-growth, the free energy of the crystal-melt interfa- 
cial free energy is an object of considerable interest across a number of scientific disciplines, especially 
in the materials-, colloid- and atmospheric sciences. Over fifty years ago, Turnbull observed that 
the interfacial free energies (scaled by the mean interfacial area per particle) of a variety of metal- 
lic elements exhibit a linear correlation with the enthalpy of fusion. This correlation provides an 
important empirical "rule-of-thumb" for estimating interfacial free energies, but lacks a compelling 
physical explanation. In this work we show that the interfacial free energies for close-packed metals 
are linearly correlated with the melting temperature, and are therefore primarily entropic in origin. 
We also show that the slope of this linear relationship can be determined with quantitative accuracy 
using a hard-sphere model, and that the correlation with the enthalpy of fusion reported by Turnbull 
follows as a consequence of the fact that the entropy of fusion for close-packed metals is relatively 
constant. 

The crystal-melt interfacial free energy, 7, defined as 
the reversible work required to form a unit area of in- 
terface between a crystal and its coexisting fluid, plays 
a central role in determining the kinetics of crystal nu- 
cleation and growth Unfortunately, direct experi- 
mental determinations of this quantity are difficult and 
exist for only a handful of materials. For most materials, 
knowledge of the interfacial free energy is obtained indi- 
rectly from measurements of crystal nucleation rates from 
undercooled fluids, from which 7 is determined within the 
approximations of classical nucleation theory. 

In a seminal 1950 paper Turnbull reported values 
of the crystal- melt interfacial free energy, 7, for a variety 
of materials, mostly metallic elements, which were ob- 
tained indirectly from nucleation rate experiments. For 
systems in which directly determined values exist with 
which to compare, the values so obtained typically are 
accurate within about 10-20%. (For example, the inter- 
facial free energy for bismuth was determined Q using 
grain boundary angles to be 61.3 xl0~ 3 J m~ 2 , as com- 
pared to the value of 54.4 xlO -3 J m~ 2 obtained from 
nucleation rate data). In order to compare the results 
for various systems, Turnbull defined a "gram-atomic" 
(or molar) interfacial free energy as the free energy of an 
interface (one atom thick) containing Avagadro's num- 
ber, Na, of atoms (or molecules): 



7 = IP 
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ture of the interface l5jj6|,|7| , but the results are quite sen- 
sitive to the nature of the assumed model. Inherent in 
these early models are two basic assumptions: 1 ) the 
solid-liquid free energy is primarily entropic in origin, 
and 2) the surface free energy is due to the increase in 
entropy associated with the enhanced structure of the liq- 
uid at the interface (i.e., the crystal is unchanged up to 
the interface). The first assumption is reasonable given 
the dominant role that packing considerations play in de- 
termining the structure and thermodynamics of simple 
liquids || and implies that a hard-sphere model should 
be adequate to describe the interfacial system. How- 
ever, the assumption 2) is at variance with molecular- 
dynamics simulations of crystal-melt interfaces of simple 
systems [p|Jl0|l that show significant structural relaxation 
(evidenced by an increase in the mean-squared displace- 
ment from the lattice sites) occurs in the crystal as the 
interface is approached. Therefore, any theory of the 
crystal-melt interfacial free energy must include a realis- 
tic description of both the solid and fluid in the interfacial 
region. 

Recently, we have determined via molecular-dynamics 
computer simulation, the structure [ fiof and solid-liquid 
interfacial free energy (n] for a system of hard-spheres. 
For this system, which freezes into a face-centered-cubic 
(fee) crystal structure, 7 was determined to be slightly 
anisotropic with an orientationally averaged value of 



The data for 7 was found to exhibit a strong correla- 
tion with the latent heat of fusion. Empirically, Turnbull 
found 



7/i -s 



0.61 kT rn /a 2 



(3) 
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where the Turnbull coefficient, Ct, was found to be ap- 
proximately 0.45 for metals (especially close-packed met- 
als) and 0.32 for many nonmetals. 

There have been attempts to explain this empirical re- 
sult through the analysis of simple models for the struc- 



where k is Boltzmann's constant, T m , is the melting tem- 
perature and a is the hard-sphere diameter. (The scaling 
with kT rn is a consequence of the fact that the phase be- 
havior in any hard-core systems is purely entropic.) This 
value is consistent with a value of 0.55±0.02 kT m /a 2 
obtained from an analysis Jl2j (using nucleation the- 
ory) of the experimental crystallization kinetics of sil- 
ica spheres, a system well described by a hard-sphere 



model. Note that, jhs is considerably lower that that 
obtained for a fluid at a structureless hard wall, which 
was recently calculated [Q at the melting density to be 
1.99±0.18 kT m /a 2 , indicating that the entropy increase 
due to the relaxation of the crystal structure near the 
interface plays an important role in determining the in- 
terfacial thermodynamics. 

For the hard-sphere system, the density of the solid 
at freezing is approximately per 3 = 1.04, independent of 
temperature [ fl3[ , therefore TurnbulPs gram-atomic inter- 
facial free energy can be written 

%,/R = 0.59(0.54)T m , (4) 

where R is the gas constant and the value in parenthesis 
is that obtained using the value of jhs determined from 
nucleation experiments [jfj). Eq. 4 suggests that if the 
interfacial free energy of fec-forming materials is well de- 
scribed by a hard-sphere model, then one should see a 
linear correlation between j/R and the melting temper- 
ature with a slope in the range 0.5 to 0.6. 
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FIG. 1. Gram-atomic interfacial free energy (scaled by the 
gas constant to give units of Kelvin) for a variety of elemental 
systems. Fcc-crystal-forming systems are shown as filled cir- 
cles and values for non-close-packed crystal formers are shown 
as open diamonds. The dashed line is a line of slope 0.50 rep- 
resenting the best linear fit passing through the origin to the 
data for the fee crystal formers. 

To test this hypothesis, we plot in Fig. [I], 7/i? as a func- 
tion of melting temperature for a variety of elemental 
materials using TurnbulPs original data. Although all of 
the data exhibits a correlation with T m , the correlation 
for the fcc-crystal forming metals is linear with a slope of 
0.50, whereas that for non-fee forming materials exhibits 



significant scatter. This value of the slope is about 20% 
below that predicted by the direct hard-sphere results, 
but less than 10% below that based on the nucleation 
rate result, which is probably a more relevant compar- 
ison, since TurnbulPs data was also so obtained. That 
the attractive forces contribute only about 10% to the 
interfacial free energy is consistent with results from sim- 
ulation H[ll| and density-functional theory [jl5j for the 
Lennard-Jones system, a prototypical fec-forming model 
potential. It should be noted that Turnbull also reported 
H a correlation of 7 with T m , but, due to the scatter in 
the overall data, it was rejected as the basis for an em- 
pirical rule in favor of the correlation with the enthalpy 
of fusion. However, he did observe that the correlation 
with T m was sensitive to the " complexity" of the crystal 
structure. This is evident in the data shown in Fig. [l], 
which shows a strong linear correlation with T m , for fec- 
forming metals, whereas that for the non-close-packed 
materials is much weaker. The hard-sphere interaction 
does not have the long-range forces necessary to mechan- 
ically stabilize a non-close packed crystal structure. As a 
consequence, the interfacial thermodynamics of systems 
that freeze into open crystal lattices will not be well de- 
scribed by a purely entropic model and a simple linear 
scaling with the melting temperature is not expected. 

For fee- (and probably hep-) forming materials, it is rel- 
atively straightforward now to understand the empirical 
correlation of the interfacial free energy with the heat of 
fusion. For the close-packed systems studied by Turnbull, 
the entropy of fusion, A f us S, is very nearly constant with 
an average value very close to that for the hard-sphere 
system, where Af us S = 9.7 J/(mole K). At the melting 
point, equilibrium requires that Af us H — T m Af us S, so 
the enthalpy of fusion should scale nearly linearly with 
the melting temperature for these systems. So if Af us H 
scales approximately linear with T m and -y^s is propor- 
tional to T m , it then follows that 7^ s will exhibit strong 
linear correlation with Af us H, and thus TurnbulPs rule 
obtains. 
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